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Abstract 



We propose a ^-linear scheme for the numerical solution of the quasi-static 
Maxwell-Landau-Lifshitz-Gilbert (MLLG) equations. Despite the strong non- 
linearity of the Landau-Lifshitz-Gilbert equation, the proposed method results 
in a linear system at each time step. We prove that as the time and space steps 
tend to zero (with no further conditions when 6 £ (i, 1]), the finite element so- 
lutions converge weakly to a weak solution of the MLLG equations. Numerical 
results are presented to show the applicability of the method. 
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1 Introduction 

The Maxwell-Landau-Lifshitz-Gilbert (MLLG) equations describe the electromag- 
netic behavior of a ferromagnetic material. In this paper, for simplicity, we suppose 
that there is a bounded cavity Del 3 (with perfectly conducting outer surface dD) 
into which a ferromagnet D is embedded. We assume further that D\D is a vacuum. 
We will consider the quasi-static case of the MLLG system. Letting := (0, T) x D 
and Dt '■— (0, T) x D, the magnetization field m : — > § 2 , where § 2 is the unit 
sphere in M 3 , and the magnetic field H : Dt — > ffi 3 satisfy 

m t = Aim x H cS — \ 2 m x (m x H cS ) in D T , (1.1) 
fioH t + <rV x (V x H) = -fi m t in D T , (1.2) 
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in which Ai 7^ 0, A2 > 0, cr > 0, and /io > are constants. Here m : Dt — » R is the 
zero extension of m onto Dt, i.e., 



m(t, a;) 




(t, x) e Dt 
{t,x)eD T \D T . 



For simplicity the effective field -ff e ff is taken to be i/ e fr — Am + H. 
The system f ll . 1 j) ( TTT2|) is supplemented with initial conditions 

m(0, .) = m in D and H(0, .) = Jf in 5, (1.3) 

and boundary conditions 

<9 n m = on (0, T) x &D and (Vxi/)xn = on (0, T) x &D. (1.4) 

The equation ( II. ip is the first dynamical model for the precessional motion of a 
magnetization, suggested by Landau and Lifshitz [12J. The existence and uniqueness 
of a local strong solution of (ll.ll) -( jl~4l) is shown by Cimrak [8]. He also proposes [7] 
a finite element method to approximate this local solution and provides error estima- 
tion. 

Gilbert introduces a different approach for description of damped precession in [9] : 

Xiirit + X 2 m x m t = fxm x H e s, (1.5) 

in which \i = Xf + X 2 ,. A proof of the equivalence between (11.51) and (II. ip can be 
found in [13]. It is easier to numerically solve ( 11. 5ft than ( 11. ip because the latter has 
a double cross term, namely m x (m x -ff e ff)- 

Instead of solving (ll.ll) - (ll.4p . Bahas, Bartels and Prohl |2] propose an implicit 
nonlinear scheme to solve problem (ll.2p - (ll.5p . and prove that the finite element 
solution converges to a weak global solution of the problem. Their method requires a 
condition on the time step k and space step h (namely k = 0(h 2 )) for the convergence 
of the nonlinear system of equations resulting from the discretization. 

Following the idea developed by Alouges and Jaison [1] for the Landau-Lifshitz- 
Gilbert (LLG) equation (11.51) . we propose a 0-lin ear finite element scheme to find a 
weak global solution to (ll.2p -( fL5l) . We prove that the numerical solutions converge 
to a weak solution of the problem with no condition imposed on time step and space 
step as 9 e (|, 1]. It is required that k = o{h 2 ) when 6 G [0, |), and k = o{h) when 
9 — 2- The implementation aspect of the algorithm is reported in (13] where no 
convergence analysis is carried out. 

The paper is organized as follows. Weak solutions of the MLLG equations are 
defined in Section [2J We also introduce in this section the ^-linear finite element 
scheme. Some technical lemmas are presented in Section [3J In Section HJ we prove 
that the finite element solutions converge to a weak solution of the MLLG equations. 
Numerical experiments are presented in the last section. 
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2 Weak solutions and finite element schemes 

Before presenting the definition of a weak solution to the MLLG equations, it is 
necessary to introduce some function spaces and to assume some conditions on the 
initial functions mo and H . 

The function spaces H 1 (D,M 3 ) and H(curl; D) are defined as follows: 

H^R 3 ) = lue L 2 (D,R 3 ) : — G L 2 (D,R 3 ) for % = 1,2,3. \ , 

H(curl; 5) = |w G L 2 (£),R 3 ) : V x u G L 2 (Z),R 3 )} . 

Here, for a domain Q C R 3 , L 2 (f2,R 3 ) is the usual space of Lebesgue squared inte- 
grable functions defined on Q and taking values in R 3 . Throughout this paper, we 
denote 

(■> -)n '■= (•> -)L2(n,R3) and II ■ \\n ■= || • ||L2(n,R3). 

In order to define a weak solution of MLLG equations, we assume that the given 
functions m Q and H Q satisfy 

m G H 1 (D, R 3 ), |m | = 1 a.e. in D and H G H(curl; D). (2.1) 

For physical reasons (see [10]), these initial fields must satisfy 

div(H + X-D m o) — in D and (H + Xd^io) • n = on (2.2) 

Since equations (11.11) and (II. 5p are equivalent (a proof of which can be found 
in [13]), instead of solving (ll.ip -( TO|) we solve fl 1 .2|) ( TT75|) . A weak solution of the 
problem is defined in the following definition. 

Definition 2.1. Let the initial data (m ,H ) satisfy (12.11) and (12. 2p . Then (m,H) 
is call a weak solution to (ll.2p - (ll.5p if, for all T > 0, there hold 

1. m G H^-D^R 3 )) and \m\ = 1 a.e. in D T ; 

2. H,H tl V x H G L 2 (5 T ,R 3 ); 

3. for all </) G C°°(D T ) and C G C°°(D T ), 

Ai (rot, 0) Dt + A 2 (m x m t , 0) Dt = \i (Vm, V(m x 0))^ + /i (m x if , 0) Dt 

(2.3) 

and 

/i (H t , C)5 T + o- (V x H , V x C,)d t = -to (™>t, C) d t , (2-4) 
where n = X\ + \\; 
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4- in the sense of traces there holds 

m(0,-) = m o , (2.5) 

5. for almost all T' G (0, T), 

8{T') + A 2 ^- 1 ||m t ||^ T/ + X 2f i- 1 \\H t \\l Ti + 2^V||V x Hf^ < 8(0), (2.6) 
where 

8{T>) = \\Vm(T>)\\l + ||iT(T')||| + A 2 ^Vo V||V x H(T')\\%. 

We next introduce the ^-linear finite element scheme which approximates a weak 
solution (m, H) defined in Definition 12. 11 

Let T/i be a regular tetrahedrization of the domain D into tetrahedra of maxi- 
mal mesh-size h, and let T^|_d be its restriction to D C D. We denote by Nh '■= 
{x\, . . . ,xn} the set of vertices and by M. h ■= {ei, . . . , eu } the set of edges. 

To discretize the LLG equation H2.3I) . we introduce the finite element space C 
e^AM 3 ) which is the space of all continuous piecewise linear functions on Th\D- A 
basis for can be chosen to be (4> n )i<n<N, where <p n (x m ) = 5 n , m - Here 5 nym stands 
for the Kronecker symbol. The interpolation operator from C°(.D,M 3 ) onto is 
denoted by Jv fe , 

N 
n=l 

To discretize Maxwell's equation (12.41) . we use the space of lowest order edge 
elements of Nedelec's first family [14J. It is known [T4] that Y^ is a subspace of 
H(curl; D) and that the set {t/> x , . . . , i^ M } is a basis for Y^ if it satisfies 

if> q G {V : D ->■ R 3 | iP\ K (x) = a K + b K xx, a K , b K G M 3 , G T h }, 

(2.7) 

Xp q -Tpds = 5 qp} 

where t p is the unit vector in the direction of edge e p . We also define the following 
interpolation operator Jy h from C°°(D) onto Y^, 

M 

I Yh (u) = J2 u ^v V W GC°°(5,M 3 ), 

q=l 

where u q = j u ■ T q ds. 

Fixing a positive integer J, we choose the time step k to be k = T/J and define 
tj = jk, j = 0, • • ■ , J. For j = 1, 2, . . . , J, the functions m(tj, •) and H(tj, ■) are 
approximated by rnf 1 G and G Yh, respectively. 



4 



We define the space 



by 



w e Y h | w(x n ) ■ m^\x n ) 



0, n = l,...,N. . 



and denote 



if 



0+1/2) 



if 



(i+i) 



ff 



(j) 



and d t H 



{3+1) 
h 



H 



Algorithm 2.1. 

Step 1: S^t j = 0. Choose = Jv h mo -H^ = h h H . 



Step 2: Fznd (v 



(i+l) rr(j+l)\ 



x Yft satisfying 



D 



r (i + l/2) 



, W h 



1 



D 



(2.8 



and 



to ( d t H^ +1 \ (l + a(vx H^ 2 \ V x C h 



D 



D 



(2.9) 



Step 3: Define 



m 



(3+1), 



JV 



m^(a; n ) + £;t> 



(i+i), 



a; 



n=1 im^^+^r 1 ^^, 

Step 4: S'et j — j + 1, and return to Step \2. 1\ if j < J. Stop if j = J. 



The parameter in (12. 8p can be chosen arbitrarily in [0, 1]. The method is explicit 
when 8 = and fully implicit when 6 = 1. 

By the Lax-Milgram Theorem, for each j > there exists a unique solution 

x Y h of equations (J23MZ3). Since ™° } ' 



f,,^' +1 ) wO+ih r- 



O'+i), 



' (^n) 



m h '\ X n l 



1 and 



for all n 



, iV , there holds 



m\i\x n ) + kv 



(3+1), 



tC 7 



> 1. 



(2.10) 



Therefore, the algorithm is well defined. There also holds 
n = !,-•■, N. 



in 



0+i) 



1 for 
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3 Some technical lemmas 

In this section we present some lemmas which will be used in the rest of the paper. 
We start by recalling the following lemma proved in [3]. 

Lemma 3.1. If there holds 

/ V0j • V0j dx < for all i, j G {1, 2, • • • , J} and i j, (3.1) 

J D 

then for all it6 V/, satisfying \u(xi) \ > 1, I = 1, 2, • • • , J, there holds 

5. (3.2) 

Condition ( 13. ip holds if all dihedral angles of the tetrahedra in Th\D are less than or 
equal to 7r/2; see [3]. In the sequel we assume that ( 13 .ip holds. 

The next lemma defines a discrete L p -norm in Yh which is equivalent to the usual 
L p -norm. 

Lemma 3.2. There exist h-independent positive constants C\ and C2 such that for 
all p G [1, 00] and u G there holds 





*"(r) 




L 


dx < j Vu 






J D 



N 



c4u\\i m <h d Y,n*n)\ p <c4u\\ p LP 



(«)' 



n=l 



w/iere tt C M d ; d= 1,2,3. 

A proof of this lemma for p = 2 and d = 2 can be found in [111 Lemma 7.3] or [61 
Lemma 1.12]. The result for general values of p and d can be obtained in the same 
manner. 

The following lemma can be proved by using the technique in [TT1 Lemma 7.3] . 

Lemma 3.3. There exists an h-independent positive constant C such that for each 
tetrahedron K G T/, and v G Yh there holds 

\\v(x)\ - \v(xi)\\ < Ch\Vv(x)\ for all x G K, 

where {a?i}i=i,2,3 o^e the vertices of K. 

Finally the following lemma is elementary; the proof of which is included for 
completeness. 



Lemma 3.4. The solutions [m^ j 



rithra \2.1\ satisfy 



0, 1, • • • , J, obtained from Algo- 



rn 



(i+i), 



m h ] \x n ) 



< 



v 



Cj+i), 



Vn = l,2,...,JV, j = 0,...,J. (3.3) 
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Proof. By using the definition of m^ +1 \ the property m^\x n ) ■ v^ +1 \x n ) = 0, and 
the identity 

\m%\x n ) + kv% +1 \x n )\ = Jl + k*\v<Z +1) (x n )\* 

we obtain 

2 



.0')/ 



A- 



A- 



m^' ) (a ;n ) + A:^ +1) « 



m^a^l - \m%\x n ) + fc V r ij K)l) + toj? +1J (*» 



.(i+i)/ 



(j+i), 



2 + 2A; 2 




2 -2 v /l + P|^' +1) (^ n )l 2 


A; 2 (l + ^l-y^K)! 2 ) 





l + fc 2 lt>g +1) (g w )l 2 -l 

A; 2 a/i + A; 2 |^ i+1) (cc n )| 2 



Using the fact that 
we deduce 

m{f +1) (x„) - m^^Xn) 



2< Vl + fc 2 rf' +1) WI 2 + l 



< 



A: 2 



proving the lemma. 



□ 



In the following section, we show that our numerical solution converges to a weak 
solution of the problem f ll.2p - fll.5p . 



4 Existence of weak solutions 

The next lemma provides a bound in the L 2 -norm for the discrete solutions. 
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Lemma 4.1. The sequence I ijn^ , , ) \ produced by Alaorithm \2.1\ 

K J 3=0,1) ••• ,J 

satisfies 



j'-i j-i 

r(*+i)iia 
Id 

i=0 i=0 

i-i 



+ C J] 112, + £ fcll^ 

i=0 

+ 2 Mo -VX:fc||Vxlff 1 / 2 )|||<^ (4.1) 



i=0 



where 
and 



4 J) = HVmW||2, + ||/1?|| 2 , + A 2 ^Vo V||V x If^H 2 -, 



c 



Aa^ 1 - (1 - 29)C 1 kh~ 2 , fl€[0,|), 
m which C\ is a positive constant which is independent with j , k and h. 

Proof. Choosing = in (ESD and C h = H^ +1/2) in (Ell), we obtain 

(4.2) 

^d t \\H^% + a\\V x Ht 1/2 % = ^frUf +1/2) ) 5 . (4-3) 

Multiplying /i/1q 1 to both sides of (14.31) and adding the resulting equation to (14. 2p . 
we deduce 



A.ll^llL + fc^llV^H 2 , + ^d t \\H ( i +1 % 



+ m -V||V x Ht l/2) \\% = -|i(Vmf,V,f + \. (4.4) 



D 



Since + kv^ +1) e Y h and 



mg + fc^' +1) 
|m£° + fctf<f +1 >| 



it follows from ( I2.10P and Lemma 13.11 that 

\\V m ty D <MrnP+kvt 1} W D - 

Equivalent ly, we have 

llVm^H 2 , < ||Vm^||U^I|V^r i) || c + 2fc(Vm?^V^' +1) } i3 . (4.5) 
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Equality (14. 4p is used to obtain from (14. 5[) the following inequality 

II D 



WVrn^Wl < ||Vm?||L - k\2Q - 1)|| Vv^Wl ~ 2k\ 2 ^\\v^^ 



-kd t \\H^\\l-2k^a\\VxHt 1/2 %- 

Hence, 



\\Vmt 1) \\l+ ||^ +1) ||| + 2a 2 ^||^ +1) ||U2^ V||V x H^ 2 % 

+ k\2Q - 1)||V^ +1) ||L < \\Vmff D + Hfr^llJ. (4.6) 

Next choosing C, h = d t H^ +1 ^ in equation (12. 9p . we obtain 

2k^\\d t H^ 1 % + a\\V x =a\\V x H«% 

The term -2kfM (v% +1) , d t H^ +1) ^_ can be estimated by 

-2kno (v { r\d t Ht 1] ) 3 < *K>||t# +1) ||2, + k» \\d t H« +1) \\%. 
Therefore, we deduce 

^ ||^ +1) ll5 + x H« +1 % < a\\V x HfWl + k^\\vt +l) f D . (4.7) 

Multiplying A2/U~ 1 / u ~ 1 to both sides of (14. 7p and adding the resulting equation to 
(14. 6p . we obtain 

||Vm? +1 >||2,+ + x H^IH 

+ kX^-^m + kX 2 ^\\d t H^\\l + 2kK V||V x ffr /2) ||| 
+ fc 2 (2^-l)||V^ +1) ||? ) 
< ||Vm^|| 2 D + ||flf ||| + V-Vo V||V x H%%. 
Replacing j by i in the above inequality and summing over i from to j — 1 yield 



PWl + ll^lll + A 2 ^VoV||V x frW||l + A 2/ ^ g ^ll^ll 2 , 

£ fcii *h r> in + 2^ v g k ii v x h m 

i=0 i=0 

(i+l)||2 
h WD 

i=0 



+ A; 2 (20 - 1) £ || VuJ 

i=0 

< || Vm° h \\l + ||if°||| + A 2/ i-Vo^l|V x fl°|||. (4.8) 
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When 9 G [~, 1], the term k 2 (28 — 1) X^=o ll^ t, h +1 ^lli» * s nonnegative. Hence, from 
inequality (14. 8 j) we obtain ( 14.1 ft where C = A2yU _1 . When G [0, ~), using the inverse 
estimate we obtain 



j'-i j'-i 

,(m)||2 ^ ;„2^ fl 1\ IIV7..( i + 1 )l|2 



^ 2 /r 2 (20 - 1) J] ||„< i+1 > II 2 , < A; 2 (20 - 1) J2 II W, 



ID) 



(4.9) 



i=0 



i=0 



where, Ci is a positive constant which is independent with j, k and h. Hence, from 
inequality (14. 8 p we obtain (14.11) where C = A 2 /i~ 1 — C\kh~ 2 {l — 29). This completes 
the proof of the lemma. □ 

Remark 4.2. The constant C in the above lemma is positive when 6 [1/2,1], When 
9 G [0, 1/2) the additional condition k = o{h 2 ) assures us that C is positive when h 
and k are sufficiently small. This condition will be required later in the following 
lemma and theorem. 



The discrete solutions m£ , and H^' constructed via Algorithm 12.11 are 
interpolated in time in the following definition. 

Definition 4.3. For each t G [0,T], let j G {0, J} be such that t G [tj,t j+ i). We 
define for t G [0, T] and x G D 



m hjk (t,x 
* 

v h ,k{t,x 
H h)k (t,x 

H hjk (t,x 
and Hj lk {t,x 



t - t 



j mi j+1) (x) + 



tj+i — t 



m^\x), 



m 



0), 



x 



V 



(i+i), 



x 



t-t 



H 



(j), 



x 



- H 



CH-i), 



x 



Hf{x) 



The following lemma shows that {m h k }, {m hk } and {v h k } converge (up to the 
extraction of subsequences) as h and k tend to 0. 

Lemma 4.4. Assume that h and k go to with a further condition k = o(h 2 ) 
when 9 G [0, ~) and no condition otherwise. There exist m G M 1 (D T ,'R 3 ) and 
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H G B 1 (0,T,h 2 (D)) such that Vxff fre/onas to L 2 (L> T ) and 

rrih,k ~+ m strongly in ~L 2 (Dt), 



dm htk 
dt 



m t weakly in ~L 2 (Dt) 



— ^ m t weakly in ~L 2 (Dt), 
m^ k — > m strongly in L?(Dt), 
\m\ = 1 a.e. in Dt, 
H h>k -»> H weakly in H 1 (0, T, L 2 (£))), 
V x ffft^ — ^ V x weakly in h 2 (D T ), 
and V x iT^ fe ^VxH weakly in h 2 (D T ). 



4.10) 
4.11) 

4.12) 
4.13) 
4.14) 
4.15) 
4.16) 
4.17) 



Proo/. 

Proof of (OOj) and ( HTTj) : 

Our goal is to prove that {m^k} is bounded in H 1 (Dr,IR 3 ) and then use the 
Banach-Alaoglu Theorem. We note from Definition 14.31 that it suffices to prove that 



n (i) n 

Vm^|| D < c, 



j-i 
i=o 



m 



(i+i) 



m 



(?) 



< c, 



(4.18) 
(4.19) 

(4.20) 



where the generic constant c is independent of j, h, and k. Indeed, it follows from 
Definition 14.31 and the Cauchy-Schwarz inequality that 



j-i 

m hik \\ 2 DT <ck^ 

j=0 

j-i 

3=0 
2 J -I 

= £* 



\m 



(j+i) 



+ 



D 



m 



U) 



D 



Vm 



O'+i) 



<9f 



(i+i) 



m 



D 



Vm 



(?) 



D 



D 



In order to prove (14. 18ft we note that for every x G D there are at most 4 basis 
functions ni , 0„ 2 , 0„ 3 and 0„ 4 being nonzero at x. This together with (x ni )\ = 1 
and E?=i0n l («) = 1 yields 



\m^(x)\ 2 



i=l 



m^'(x r 



x 



< 1. 



(4.21) 
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This implies (I4.18p with a constant c = l-D) 1 / 2 where \D\ is the measure of the 
domain D. 

Inequality (14. 19ft is proved in Lemma 14. 11 In order to prove inequality (14.201) , we 
note that Lemma 13.41 and Lemma 13.21 imply 



— m yi 



< c 



v 



(i+i) 



D 



D 



By using this inequality, Lemma 14.11 and Remark 14.21 we deduce 



j-i 

3=0 



m 



(i+i) 



m 



0") 



A- 



D 



J-l 

<kj2 c \H 

3=0 



(3+1) 



< C. 



V 



The Banach-Alaoglu Theorem implies the existence of a subsequence of {mh^k} 
which converges weakly to a function m £ ET^-Dy) as k and h tend to zero. This 
implies (T4~T0D and (OTjl . 

Proof of (ESI): 

From (14. ip and and Remark 14. 2[ it is straightforward to show that {vh,k} is 
bounded in ~h 2 (D T ). Hence, there exists a subsequence of {v hyk } which converges 
weakly to a function v £ L 2 (Z>r). The problem reduces to proving that m t equals 
v in L 2 (Dt). In order to show this we choose for each ij> £ h 2 (Dj<) a sequence 
{V'i} £ OTl^r) converging to i/j in L 2 (Z>r) as i tends to infinity. We then have 



(m t - v, if)) 



< \ {m t -v,^ l -^) DT 



dm 



Of 



D T 



+ 



dt 



Vh,k,lf>i 



D T 



+ | (v h ,k - V, lf)i) 



D T 



< \\mt-v\\ DT \\i/>i-ip\\ DT 
dm hjk 



dm hyk 



+ 



- Vh k 



ot 

: Ti + • • ■ + T 4 . 



L!(D T ) 



illL°°(D T 



(4.22) 



By letting h, k — > and then z — >■ oo we have Tj — > for z = 1, 2 and 4. It remains to 
show that T 3 — y 0. 



It is clear from the definition of mjf in Algorithm 12.11 that 



m^ +l \x n )-m^{x n )-kv^>{xn) = m%> (x n ) + kv%^> (x n 



.Cj+i) 



.0')/ 



- 1. 



(4.23) 
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It easily follows from 



1 and v { i +1 \x n ) ■ m%>(x n ) = that 



(J), 



m^(x n ) + kv^\x n 



(j+i). 



The above inequality and ( 14.231) yield 



~ 2 



v 



(j+i). 



+ 1. 



rajf +1) (a^) - mj?'(a: n ) 



A- 



(3+1), 



- 2 



(j+i), 



By using Lemma 13.21 we deduce 
-(*) - Wfc,fc(*) 



dm htk 



Ot 



< ck \\v 



L 1 (D) 



for £ G [tj, tj 



Integrating both sides of this inequality with respect to t over an interval [tj,tj + i) 
and summing over j from to J — 1 yield, noting the boundedness of {||vh,A;||_D r }, 



drrih k 



dt 



- Vh k 



< ck \\vh,k\\ D < ck — > as ft, k — >■ 0. 



Li(D T ) 



Thus T 3 -> as ft,, fc -» and z -)> oo. It follows from (02"]) that 

|(m t -^) DT |=0 V^GL 2 ^). 

This proves (14712]) . 
Proof of ([USD : 

It is clear from the definition of ra hk and m^fc that for t G [i^ij+i) there holds 



\m h , k (t)-m hk (t)\\ D 



(t-tj 



m 



(i+i) 



(j) 



k 



< k 



D 



d(m h , k (t,x)) 



dt 



D 



Integrating both sides of this inequality with respect to t over an interval [tj,tj + \] 
and summing over j from to (J — 1) yield 



\ m h,k-rn Kk \\ DT < k 



dm 



h,k 



dt 



< ck ->■ as ft, k -> 0. 



The above result and (IPO]) imply fTCTg]) 
Proof of (TCMj) : 



Using Lemma [3731 and noting that \m^\x n )\ = 1 for n = 1, • • • , N, we deduce 



ImjftaOI-l 



<Cft 2 



Vmf(a;; 



for all a; G D. 
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Integrating both sides of the above inequality on [tj,tj+x) x D, using Lemma \A. II and 
noting Remark I4.2[ we obtain 



1 — \m^\x 



dx dt < ch 2 



Vm 



< <M 2 



Hence 



D T 



We infer from (14. 13ft that 



1 — \m h k \ | dx dt — > as /i, — > 0. 



Iml = 1 a.e. in Z>r- 



Proof of (T4TT5D. (OS]) and (T4TT7j): 

By using the same arguments as above, we obtain these results, completing the 
proof of the lemma. □ 

We are now able to prove the main result of this paper. 

Theorem 4.5. Assume that h and k go to with the following conditions 



'k = o(h 2 ) when < 6 < 1/2, 

k = o(h) when 9 = 1/2, 

no condition when 1/2 < 6 < 1. 



(4.24) 



Then the limits (m, H) given by Lemma 4-4 ^ s a weak solution of the MLLG equa- 
tions 0-(jZ3D- 

Proof. For any G C°°(D T ), C G C°°(D T ), and t G [tj,t j+1 ), we define 
w htk (t, •) := h h (rn^ h x 0(t, •)) and Cfc(*, : = h h (C(t, •))■ 

In equations (12. 8ft and (12. 9ft . replacing iu^ and C h by Wh, k (t) and C/i(*)> respectively, 
and using Definition I4.3[ we rewrite (I2.8l) - (l2.9p as 

- Ai (m Kk (t) x v h>k (t), w h>k (t)) D + A 2 (v h)k (t), w h>k (t)) D 

= -ix (V(mJ~ fc (t) + k6v h)k (t)), Vw h)k (t)) D + n (H h>k (t), w h>k (t) 



D 



and 



IMi { C h {t) ) + a { V x ff M (t), V x C h (t) )_. = -//o Ch(*))5 
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Integrating both sides of these equations with respect to t over an interval [tj,tj+i) 
and summing over j from to J — 1 yield 



^1 ( m h,k X V h,k, w h,k) Dr + A 2 (V hjk , Wh,k) DT 

= -fi <V(m^ fc + k6v htk ), Vw h<k ) DT + fi (H h}k , w K 



and 



cx(V x H h>k ,V x C, 



h ~ 



"MO \Vh,k, ChlDn 



(4.25) 



(4.26) 



In order to prove that m and H satisfy (12. 3p and (12.41) . respectively, we prove 
that as h and k tend to there hold 



(mj~ fc x « fc>fc , w^) (mx m t ,mx <f>) Dt , 
(« hifc , w Kk ) Dt ->■ (m,, mx (j)) DT , 
(Vrn htk ,Vw htk ) DT -)• (Vm,V(m x 0)) Dt , 
fc (^v h;k ,Vw h!k } DT -»■ 0, 



and 



\ <9t 



On 



Vxif M) Vx( ft \ ->(Vx JT,VxC)n 



(Vh,kXh)5 T -> ( m tX)D T - 

We now prove (I4.27P and ()4.3Up ; the others can be obtained in the same manner. 
Using the triangular inequality and Holder's inequality, we estimate 



(4.27) 
(4.28) 
(4.29) 
(4.30) 
(4.31) 

(4.32) 

(4.33) 
(4.34) 



lh,k 



m ft,fc x v Kk, Wh,k) D - (mx m t , mx (j)) 



On 



as follows 

Ih,h < 



[m h k x v hjk , w hjk - m h k x cj>) 
{m~ k - m) x v hjk , m x cf>) 



Dn 



+ 



( m h,k x v h,k, {m htk - m) x 4>) Dt 
m x (v hjk - m t ), m x <f>) DT 



h,k / ' ' ")«' ' ' ~r i q 

< \\rn Kk \\^ {DT) \\v h)k \\ DT \\w h)k - m hk x (f>\\ DT 

+ ll m ft,fcllL°°(D T )||f ft,fe||z) T ||wi - m h,k\\D T \\<t>\\h°°(D T ) 
+ ||m - m^ Jfe ||£» T ||«/ l)fc ||£) T ||0||L«>(D T ) 

+ \\Vh t k - m t \\D T \\<i>\\L°°(v T ) 

< c (\\w h> k - m Kk x 0|| Dt + ||m - m^ k \\ DT + ||« ft)A . - m t ||£, T ) 
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where we have used (I4.2ip and Lemma I4.1[ noting Remark 14.21 The interpolation 
operators Iy h and Iy h have the following properties (see e.g., [5] and [14"] ) 

(4.35) 



\\ m h,k X ^ ~ W h,k\\t, 2 ([0,T],B 1 (D)) — C^ll m / l ,A:llEI 1 ( J D T )||0||w 2 .°°(D T ) ; 

IIC(*j - < h (t)\\n+ liv x (C(t) - C fc W)lb ^ ^llv 2 CII 5 . 



This implies 



lim 4 fe = 0, 



proving f|4.27j) . 

In order to prove (14.301) we first note that 

\\Vw h;k \\ DT < ||V(m h>fc x - w h>k )\\ DT + || V(mr (fc x 0)|| Dt 

< c/i ll m Mll H1 (^T)ll < / > llw2.-(D T ) + l|Vm^| Dr ||V0|| L °°(r> T ) 

< C|| 011^2,00 (£> T ), 

where we have used (14.351) and the boundedness of H?"-^ Jln^Or)- Now using Holder's 
inequality we obtain 

k (Vv hjk , Vw Kk ) DT < c^UVv^fcH^. (4.36) 

It is straightforward from (14. 8 p that {|| Vv/^IIdt} is bounded when 9 G (|, 1]. Hence, 
taking the limit as k and h tend to in (14.361) yields (14.301) for these values of 9. 
When 9 G [0, |], using the inverse estimate we obtain 



3-1 3-i 

)||2 
Di 

i=0 i=0 



or equivalently, 

HVv^fcHoy < chr l \\v Kk \\ DT . 
Hence under the assumption (14.241) . the inequality (I4.36P becomes 

k(Sv h)k ,Vw h)k ) DT < ckh~ 1 \\v h>k \\ DT < ckh" 1 

when 9 G [0, 1/2]. Therefore, under the assumption (I4.24p there holds 

k (Vv Kk , Vw h)k } DT ->• 0. 

We now prove (12 .5p . Since m° h = Iy h (mo), the sequence {m^} converges to mo 
in ~h 2 (D) as h tends to 0. Using the weak continuity of the trace operator we obtain 
that m(0, •) = m in the sense of traces. 

Finally, applying weak lower semicontinuity of norms in inequality ( 14. ip we obtain 
the energy inequality ( 12. 6p . which completes the proof. □ 
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5 Numerical experiments 



In order to carry out physically relevant experiments, the initial fields mo, Hq must 
satisfy condition ( 12.2ft . This can be achieved by taking 

H = H* - XDm , 

where div Hq = in D. In our experiment, for simplicity, we choose Hq to be a 
constant. We solve an academic example with D — D — (0, l) 3 and 



where x = (xi,x 2 , x 3 ), x* = (x\— 0.5, x 2 — 0.5, 0) and A = (1— 2|a;*|) 4 /4. The constant 
H s represents the strength of H in the ^-direction. We compute the experiments 
for H s = 0, ±30, ±100 and ±1000. We set the values for the other parameters in (11.11) 
and (11.21) as Ai = A 2 = fio = °" = 1- 

The domain D is partitioned into uniform cubes with the mesh size h = 1/2 3 , 
where each cube consists of six tetrahedra. We choose the time step k = 10~ 3 and 
the parameter 9 in Algorithm 12. II to be 0.7. The construction of the basis functions 
for Wj^ and in this algorithm is discussed in [13]. At each iteration we need to 
solve a linear system of size (2A^ + M) x (2A^ + M), recalling that is the number 
of vertices and M is the number of edges in the triangulation. The code is written in 
Fortran90. 

The evolution of HVm^Hu, H-H^H^ and || V x H hjk \\^ are depicted in Figures [fl 
[2] and [3X respectively. Figure 0] shows that the solution satisfies condition (I2.6P in 
Definition O 

Remark 5.1. By the time this paper was written up, we learnt that Banas, Page and 
Praetorius %4}j independently solved a similar problem. They also used a linear scheme 
similar to our scheme, even though their variational formulation was different. 
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Figure 1: Plot of t H- \\Vm h , k {t)\\ D . 
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Figure 2: Plot of t H> \\H hjk (t)\\ s . 
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Figure 3: Plot of 1 1-> ||V x H 
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